Mitochondrial DNA variation in Aedes aegypti (Diptera: Culicidae) mosquitoes from Jeddah, Saudi Arabia.

Wolbachia (Hertig 1936) (Rickettsiales: Ehrlichiaceae) has emerged as a valuable biocontrol tool in the fight against dengue by suppressing the transmission of the virus through mosquitoes. Monitoring the dynamics of Wolbachia is crucial for evaluating the effectiveness of release programs. Mitochondrial (mtDNA) markers serve as important tools for molecular tracking of infected mitochondrial backgrounds over time but require an understanding of the variation in release sites. In this study, we investigated the mitochondrial lineages of Aedes aegypti (Linnaeus 1762) in Jeddah, Saudi Arabia, which is a prospective release site for the "wAlbBQ" Wolbachia-infected strain of this mosquito species. We employed a combination of comprehensive mitogenomic analysis (including all protein-coding genes) and mtDNA marker analysis (cox1 and nad5) using data collected from Jeddah. We combined our mitogenome and mtDNA marker data with those from previous studies to place mitochondrial variation in Saudi Arabia into a broader global context. Our findings revealed the presence of 4 subclades that can be broadly categorized into 2 major mitochondrial lineages. Ae. aegypti mosquitoes from Jeddah belonged to both major lineages. Whilst mitogenomic data offered a higher resolution for distinguishing Jeddah mosquitoes from the wAlbBQ strain, the combination of cox1 and nad5 mtDNA markers alone proved to be sufficient. This study provides the first important characterization of Ae. aegypti mitochondrial lineages in Saudi Arabia and offers essential baseline information for planning future molecular monitoring efforts during the release of Wolbachia-infected mosquitoes.

Adaptive evolution to the natural and anthropogenic environment in a global invasive crop pest, the cotton bollworm.

The cotton bollworm, Helicoverpa armigera, is set to become the most economically devastating crop pest in the world, threatening food security and biosafety as its range expands across the globe. Key to understanding the eco-evolutionary dynamics of H. armigera, and thus its management, is an understanding of population connectivity and the adaptations that allow the pest to establish in unique environments. We assembled a chromosome-scale reference genome and re-sequenced 503 individuals spanning the species range to delineate global patterns of connectivity, uncovering a previously cryptic population structure. Using a genome-wide association study (GWAS) and cell line expression of major effect loci, we show that adaptive changes in a temperature- and light-sensitive developmental pathway enable facultative diapause and that adaptation of trehalose synthesis and transport underlies cold tolerance in extreme environments. Incorporating extensive pesticide resistance monitoring, we also characterize a suite of novel pesticide and Bt resistance alleles under selection in East China. These findings offer avenues for more effective management strategies and provide insight into how insects adapt to variable climatic conditions and newly colonized environments.

Developing Wolbachia-based disease interventions for an extreme environment.

Aedes aegypti mosquitoes carrying self-spreading, virus-blocking Wolbachia bacteria are being deployed to suppress dengue transmission. However, there are challenges in applying this technology in extreme environments. We introduced two Wolbachia strains into Ae. aegypti from Saudi Arabia for a release program in the hot coastal city of Jeddah. Wolbachia reduced infection and dissemination of dengue virus (DENV2) in Saudi Arabian mosquitoes and showed complete maternal transmission and cytoplasmic incompatibility. Wolbachia reduced egg hatch under a range of environmental conditions, with the Wolbachia strains showing differential thermal stability. Wolbachia effects were similar across mosquito genetic backgrounds but we found evidence of local adaptation, with Saudi Arabian mosquitoes having lower egg viability but higher adult desiccation tolerance than Australian mosquitoes. Genetic background effects will influence Wolbachia invasion dynamics, reinforcing the need to use local genotypes for mosquito release programs, particularly in extreme environments like Jeddah. Our comprehensive characterization of Wolbachia strains provides a foundation for Wolbachia-based disease interventions in harsh climates.

Close Kin Dyads Indicate Intergenerational Dispersal and Barriers.

AbstractThe movement of individuals through continuous space is typically constrained by dispersal ability and dispersal barriers. A range of approaches have been developed to investigate these. Kindisperse is a new approach that infers recent intergenerational dispersal (σ) from close kin dyads and appears particularly useful for investigating taxa that are difficult to observe individually. This study, focusing on the mosquito Aedes aegypti, shows how the same close kin data can also be used for barrier detection. We empirically demonstrate this new extension of the method using genome-wide sequence data from 266 Ae. aegypti. First, we use the spatial distribution of full-sib dyads collected within one generation to infer past movements of ovipositing female mosquitoes. These dyads indicated the relative barrier strengths of two roads and performed favorably against alternative genetic methods for detecting barriers. We then use Kindisperse to quantify recent intergenerational dispersal (σ=81.5-197.1 m generation-1/2) from the difference in variance between the sib and the first cousin spatial distributions and, from this, estimate effective population density (ρ=833-4,864 km-2). Dispersal estimates showed general agreement with those from mark-release-recapture studies. Barriers, σ, ρ, and neighborhood size (331-526) can inform forthcoming releases of dengue-suppressing Wolbachia bacteria into this mosquito population.

Spatio-Temporal Modelling Informing Wolbachia Replacement Releases in a Low Rainfall Climate.

Releases of Aedes aegypti carrying Wolbachia bacteria are known to suppress arbovirus transmission and reduce the incidence of vector-borne diseases. In planning for Wolbachia releases in the arid environment of Jeddah, Saudi Arabia, we collected entomological data with ovitraps across a 7-month period in four locations. Herein, we show that mosquito presence in basements does not differ from that of non-basement areas of buildings. In modelling mosquito presence across the study sites, we found the spatial structure to be statistically significant in one of the four sites, while a significant spatial structure was found for egg production data across three of the four sites. The length scales of the spatial covariance functions fitted to the egg production data ranged from 143 m to 574 m, indicating that high productivity regions can be extensive in size. Rank-correlation analyses indicated that mosquito presence tended to persist from the dry to wet season, but that egg production ranks at locations could reverse. The data suggest that, in Jeddah, the quality of the local environment for breeding can vary over time. The data support the feasibility of dry season releases but with release numbers needing to be flexible depending on local rates of invasion.

Unveiling biogeographical patterns in the worldwide distributed Ceratitis capitata (medfly) using population genomics and microbiome composition.

Invasive species are among the most important, growing threats to food security and agricultural systems. The Mediterranean medfly, Ceratitis capitata, is one of the most damaging representatives of a group of rapidly expanding species in the family Tephritidae, due to their wide host range and high invasiveness potential. Here, we used restriction site-associated DNA sequencing (RADseq) to investigate the population genomic structure and phylogeographical history of medflies collected from six sampling sites, including Africa (South Africa), the Mediterranean (Spain, Greece), Latin America (Guatemala, Brazil) and Australia. A total of 1907 single nucleotide polymorphisms (SNPs) were used to identify two genetic clusters separating native and introduced ranges, consistent with previous findings. In the introduced range, all individuals were assigned to one genetic cluster except for those in Brazil, which showed introgression of an additional genetic cluster that also appeared in South Africa, and which could not be previously identified using microsatellite markers. Moreover, we assessed the microbial composition variations in medfly populations from selected sampling sites using amplicon sequencing of the 16S ribosomal RNA (V4 region). Microbiome composition and structure were highly similar across geographical regions and host plants, and only the Brazilian specimens showed increased diversity levels and a unique composition of its microbiome compared to other sampling sites. The unique SNP patterns and microbiome features in the Brazilian specimens could point to a direct migration route from Africa with subsequent adaptation of the microbiota to the specific conditions present in Brazil. These findings significantly improve our understanding of the evolutionary history of the global medfly invasions and their adaptation to newly colonized environments.

Global population genomic signature of Spodoptera frugiperda (fall armyworm) supports complex introduction events across the Old World.

Native to the Americas, the invasive Spodoptera frugiperda (fall armyworm; FAW) was reported in West Africa in 2016, followed by its chronological detection across the Old World and the hypothesis of an eastward Asia expansion. We explored population genomic signatures of American and Old World FAW and identified 12 maternal mitochondrial DNA genome lineages across the invasive range. 870 high-quality nuclear single nucleotide polymorphic DNA markers identified five distinct New World population clusters, broadly reflecting FAW native geographical ranges and the absence of host-plant preferences. We identified unique admixed Old World populations, and admixed and non-admixed Asian FAW individuals, all of which suggested multiple introductions underpinning the pest's global spread. Directional gene flow from the East into eastern Africa was also detected, in contrast to the west-to-east spread hypothesis. Our study demonstrated the potential of population genomic approaches via international partnership to address global emerging pest threats and biosecurity challenges.

Genomic insights into a population of introduced European rabbits Oryctolagus cuniculus in Australia and the development of genetic resistance to rabbit hemorrhagic disease virus.

The European rabbit (Oryctolagus cuniculus) is one of the most devastating invasive species in Australia. Since the 1950s, myxoma virus (MYXV) and rabbit haemorrhagic disease virus (RHDV) have been used to manage overabundant rabbit populations. Resistance to MYXV was observed within a few years of the release. More recently, resistance to lethal RHDV infection has also been reported, undermining the efficiency of landscape-scale rabbit control. Previous studies suggest that genetic resistance to lethal RHDV infection may differ locally between populations, yet the mechanisms of genetic resistance remain poorly understood. Here, we used genotyping by sequencing (GBS) data representing a reduced representation of the genome, to investigate Australian rabbit populations. Our aims were to understand the relationship between populations and identify possible genomic signatures of selection for RHDV resistance. One population we investigated had previously been reported to show levels of resistance to lethal RHDV infection. This population was compared to three other populations with lower or no previously reported RHDV resistance. We identified a set of novel candidate genes that could be involved in host-pathogen interactions such as virus binding and infection processes. These genes did not overlap with previous studies on RHDV resistance carried out in different rabbit populations, suggesting that multiple mechanisms are feasible. These findings provide useful insights into the different potential mechanisms of genetic resistance to RHDV virus which will inform future functional studies in this area.

Voltage-sensitive sodium channel (Vssc) mutations associated with pyrethroid insecticide resistance in Aedes aegypti (L.) from two districts of Jeddah, Kingdom of Saudi Arabia: baseline information for a Wolbachia release program.

BACKGROUND: Dengue suppression often relies on control of the mosquito vector, Aedes aegypti, through applications of insecticides of which the pyrethroid group has played a dominant role. Insecticide resistance is prevalent in Ae. aegypti around the world, and the resulting reduction of insecticide efficacy is likely to exacerbate the impact of dengue. Dengue has been a public health problem in Saudi Arabia, particularly in Jeddah, since its discovery there in the 1990s, and insecticide use for vector control is widespread throughout the city. An alternative approach to insecticide use, based on blocking dengue transmission in mosquitoes by the endosymbiont Wolbachia, is being trialed in Jeddah following the success of this approach in Australia and Malaysia. Knowledge of insecticide resistance status of mosquito populations in Jeddah is a prerequisite for establishing a Wolbachia-based dengue control program as releases of Wolbachia mosquitoes succeed when resistance status of the release population is similar to that of the wild population. METHODS: WHO resistance bioassays of mosquitoes with deltamethrin, permethrin and DDT were used in conjunction with TaqMan® SNP Genotyping Assays to characterize mutation profiles of Ae. aegypti. RESULTS: Screening of the voltage-sensitive sodium channel (Vssc), the pyrethroid target site, revealed mutations at codons 989, 1016 and 1534 in Ae. aegypti from two districts of Jeddah. The triple mutant homozygote (1016G/1534C/989P) was confirmed from Al Safa and Al Rawabi. Bioassays with pyrethroids (Type I and II) and DDT showed that mosquitoes were resistant to each of these compounds based on WHO definitions. An association between Vssc mutations and resistance was established for the Type II pyrethroid, deltamethrin, with one genotype (989P/1016G/1534F) conferring a survival advantage over two others (989S/1016V/1534C and the triple heterozygote). An indication of synergism of Type I pyrethroid activity with piperonyl butoxide suggests that detoxification by cytochrome P450s accounts for some of the pyrethroid resistance response in Ae. aegypti populations from Jeddah. CONCLUSIONS: The results provide a baseline for monitoring and management of resistance as well as knowledge of Vssc genotype frequencies required in Wolbachia release populations to ensure homogeneity with the target field population. Vssc mutation haplotypes observed show some similarity with those from Ae. aegypti in southeast Asia and the Indo-Pacific, but the presence of the triple mutant haplotype in three genotypes indicates that the species in this region may have a unique population history.

Adaptive Introgression across Semipermeable Species Boundaries between Local Helicoverpa zea and Invasive Helicoverpa armigera Moths.

Hybridization between invasive and native species has raised global concern, given the dramatic increase in species range shifts and pest outbreaks due to anthropogenic dispersal. Nevertheless, secondary contact between sister lineages of local and invasive species provides a natural laboratory to understand the factors that determine introgression and the maintenance or loss of species barriers. Here, we characterize the early evolutionary outcomes following secondary contact between invasive Helicoverpa armigera and native H. zea in Brazil. We carried out whole-genome resequencing of Helicoverpa moths from Brazil in two temporal samples: during the outbreak of H. armigera in 2013 and 2017. There is evidence for a burst of hybridization and widespread introgression from local H. zea into invasive H. armigera coinciding with H. armigera expansion in 2013. However, in H. armigera, the admixture proportion and the length of introgressed blocks were significantly reduced between 2013 and 2017, suggesting selection against admixture. In contrast to the genome-wide pattern, there was striking evidence for adaptive introgression of a single region from the invasive H. armigera into local H. zea, including an insecticide resistance allele that increased in frequency over time. In summary, despite extensive gene flow after secondary contact, the species boundaries are largely maintained except for the single introgressed region containing the insecticide-resistant locus. We document the worst-case scenario for an invasive species, in which there are now two pest species instead of one, and the native species has acquired resistance to pyrethroid insecticides through introgression.

Population genetics and migration pathways of the Mediterranean fruit fly Ceratitis capitata inferred with coalescent methods.

BACKGROUND: Invasive species are a growing threat to food biosecurity and cause significant economic losses in agricultural systems. Despite their damaging effect, they are attractive models for the study of evolution and adaptation in newly colonised environments. The Mediterranean fruit fly, Ceratitis capitata, as a member of the family Tephritidae, is one of the most studied invasive species feeding on many fruit crops in the tropics and subtropics worldwide. This study aims to determine the global macrogeographic population structure of Ceratitis capitata and reconstruct its potential migration routes. METHOD: A partial mitochondrial cytochrome oxidase I gene from >400 individual medflies and 14 populations from four continents was sequenced and subjected to Bayesian demographic modelling. RESULTS: The Afrotropical populations (Kenya, South Africa and Ghana) harbour the majority of haplotypes detected, which also are highly divergent, in accordance with the presumed ancestral range of medflies in Sub-Saharan Africa. All other populations in the presumed non-native areas were dominated by a single haplotype also present in South Africa, in addition to a few, closely related haplotypes unique to a single local population or regional set, but missing from Africa. Bayesian coalescence methods revealed recent migration pathways from Africa to all continents, in addition to limited bidirectional migration among many local and intercontinental routes. CONCLUSION: The detailed investigation of the recent migration history highlights the interconnectedness of affected crop production regions worldwide and pinpoints the routes and potential source areas requiring more specific quarantine measures.

Standardized molecular diagnostic tool for the identification of cryptic species within the Bemisia tabaci complex.

BACKGROUND: The whitefly Bemisia tabaci complex harbours over 40 cryptic species that have been placed in 11 phylogenetically distinct clades based on the molecular characterization of partial mitochondrial DNA COI (mtCOI) gene region. Four cryptic species are currently within the invasive clade, i.e. MED, MEAM1, MEAM2 and IO. Correct identification of these species is a critical step towards implementing reliable measures for plant biosecurity and border protection; however, no standardized B. tabaci-specific primers are currently available which has caused inconsistencies in the species identification processes. RESULTS: We report three sets of polymerase chain reaction (PCR) primers developed to amplify the mtCOI region which can be used for genotyping MED, MEAM1 and IO species, and tested these primers on 91 MED, 35 MEAM1 and five IO individuals. PCR and sequencing of amplicons identified a total of 21, six and one haplotypes in MED, MEAM1 and IO respectively, of which six haplotypes were new to the B. tabaci database. CONCLUSION: These primer pairs enabled standardization and robust molecular species identification via mtCOI screening of the targeted invasive cryptic species and will improve quarantine decisions. Use of this diagnostic tool could be extended to other species within the complex. © 2017 Society of Chemical Industry.

The Trouble with MEAM2: Implications of Pseudogenes on Species Delimitation in the Globally Invasive Bemisia tabaci (Hemiptera: Aleyrodidae) Cryptic Species Complex.

Molecular species identification using suboptimal PCR primers can over-estimate species diversity due to coamplification of nuclear mitochondrial (NUMT) DNA/pseudogenes. For the agriculturally important whitefly Bemisia tabaci cryptic pest species complex, species identification depends primarily on characterization of the mitochondrial DNA cytochrome oxidase I (mtDNA COI) gene. The lack of robust PCR primers for the mtDNA COI gene can undermine correct species identification which in turn compromises management strategies. This problem is identified in the B. tabaci Africa/Middle East/Asia Minor clade which comprises the globally invasive Mediterranean (MED) and Middle East Asia Minor I (MEAM1) species, Middle East Asia Minor 2 (MEAM2), and the Indian Ocean (IO) species. Initially identified from the Indian Ocean island of Réunion, MEAM2 has since been reported from Japan, Peru, Turkey and Iraq. We identified MEAM2 individuals from a Peruvian population via Sanger sequencing of the mtDNA COI gene. In attempting to characterize the MEAM2 mitogenome, we instead characterized mitogenomes of MEAM1. We also report on the mitogenomes of MED, AUS, and IO thereby increasing genomic resources for members of this complex. Gene synteny (i.e., same gene composition and orientation) was observed with published B. tabaci cryptic species mitogenomes. Pseudogene fragments matching MEAM2 partial mtDNA COI gene exhibited low frequency single nucleotide polymorphisms that matched low copy number DNA fragments (<3%) of MEAM1 genomes, whereas presence of internal stop codons, loss of expected stop codons and poor primer annealing sites, all suggested MEAM2 as a pseudogene artifact and so not a real species.

Identification and preliminary characterization of chemosensory perception-associated proteins in the melon fly Bactrocera cucurbitae using RNA-seq.

An investigation into proteins involved in chemosensory perception in the melon fly, Bactrocera cucurbitae (Diptera: Tephritidae) is described here using a newly generated transcriptome dataset. The melon fly is a major agricultural pest, widely distributed in the Asia-Pacific region and some parts of Africa. For this study, a transcriptome dataset was generated using RNA extracted from 4-day-old adult specimens of the melon fly. The dataset was assembled and annotated via Gene Ontology (GO) analysis. Based on this and similarity searches to data from other species, a number of protein sequences putatively involved in chemosensory reception were identified and characterized in the melon fly. This included the highly conserved "Orco" along with a number of other less conserved odorant binding protein sequences. In addition, several sequences representing putative ionotropic and gustatory receptors were also identified. This study provides a foundation for future functional studies of chemosensory proteins in the melon fly and for making more detailed comparisons to other species. In the long term, this will ultimately help in the development of improved tools for pest management.

The complete mitochondrial DNA genome of a Chloridea (Heliothis) subflexa (Lepidoptera: Noctuidae) morpho-species.

We investigate the complete mitogenome of a pheromone-trapped morpho-species of Chloridea subflexa from Brazil (initially identified by the Sanger sequencing of partial mtCOI gene) as 15 323 bp (KT598688) via next generation sequencing platform. The mitogenome has an A/T rich base composition (A: 40.4%; T: 40.3%; C: 11.5%; G: 7.8%), and included 13 protein-coding genes (PCGs), 22 tRNAs, 2 ribosomal RNAs and a putative replication region (ca. 323 bp). All PCGs start with a methionine (M) amino acid except the COI gene which has an arginine (R). The trnL2 and trn-Lys genes were partially embedded within the COII gene, while the trn-His gene was completely embedded within the ND4 gene. All PCGs ends with the "TAA" stop codon except ND3 which has a "TAG" stop codon.

Detection of the Acetylcholinesterase Insecticide Resistance Mutation (G328A) in Natural Populations of Ceratitis capitata.

Wild Mediterranean fruit fly specimens collected from various regions worldwide were screened for the glycine to alanine (Gly->Ala) point mutation (G328A) in the acetylcholinesterase enzyme, presumably causing resistance to organophosphates. We found that the single nucleotide polymorphism (SNP) responsible for this amino acid change is located at the beginning of exon 6 of the Ccace2 gene. The identification of the exact location of the SNP permitted PCR primer design around this site and direct sequencing of the corresponding genomic region. We detected the resistance allele in natural Mediterranean fruit fly populations from Brazil and Spain, but not from other sites in four continents. The known treatment history of sites suggests that the resistance buildup is linked to organophosphate application in the field. The PCR-based detection provides a screening method useful for monitoring Mediterranean fruit fly insecticide resistance in local populations and improving pest management strategies accordingly.

Molecular taxonomic identification of Dacus and Ceratitis species from Sub-Saharan Africa using mitochondrial haplotypes.

Several different taxa within the genera Dacus and Ceratitis (Diptera: Tephritidae) are important agricultural pests in Sub-Saharan Africa. Although the status of many of these taxa as distinct species and their phylogenetic relationships is unclear, it is clear that these pests use a wide range of host plants and are highly invasive. The great potential for economic damage inflicted by these pests requires the ability to make accurate and reliable taxonomic identification of specimens. However, many limitations and uncertainties are encountered when these species are examined using traditional approaches based on morphological identification techniques. We describe here the amplification and analysis of DNA sequences from the mitochondrial cytochrome oxidase II (COII)-tRNA(lys)-tRNA(ASP) genes from individuals of various Dacus and Ceratitis species and populations from Sub-Saharan Africa. The variation detected in the DNA sequences of these individuals is used both for clarification of their taxonomic status and the analysis of phylogenetic relationships of these taxa.